Steady state work fluctuations of a dragged particle under external and thermal noise 
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We consider a particle, confined to a moving harmonic potential, under the influence of friction 
and external asymmetric Poissonian shot noise (PSN). We study the fluctuations of the work done 
to maintain this system in a nonequilibrium steady state. PSN generalizes the usual Gaussian noise 
and can be considered to be a paradigm of external noise, where fluctuation and dissipation originate 
from physically independent mechanisms. We consider two scenarios: (i) the noise is given purely 
by PSN and (ii) in addition to PSN the particle is subject to white Gaussian noise. In both cases we 
derive exact expressions for the large deviation form of the work distribution, which are characterized 
by the time scales of the system. We show that the usual steady state fluctuation theorem is violated 
in our model and that in a certain parameter regime large negative work fluctuations are more likely 
to occur than the corresponding positive ones, though the average work is always positive. 
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I. INTRODUCTION 

An understanding of the fluctuation properties of sys- 
tems away from thermal equilibrium plays an important 
role for a statistical mechanical theory of these systems. 
Nonequilibrium fluctuations can have significant effects 
in small scale systems from physics to biology, which be- 
come more and more important in technological inno- 
vations. However, from a fundamental theoretical point 
of view, not much is known about the general aspects of 
nonequilibrium fluctuations, contrary to our understand- 
ing of equilibrium fluctuations. 

The simplest generalization of the equilibrium state 
is a nonequilibrium steady state (NESS), which arises 
physically due to a balance between driving forces and 
dissipative forces acting on the system. One of the sim- 
plest nonequilibrium systems that can be maintained in 
a NESS is a Brownian particle, e.g., a spherical colloidal 
particle in a fluid, confined to a harmonic potential which 
moves with constant velocity. In addition to the poten- 
tial force the particle is subject to friction and thermal 
noise from the surrounding fluid. In this system the driv- 
ing is due to the time-dependent force from the moving 
potential and the dissipation is due to the friction of the 
particle in the fluid. In order to maintain the NESS one 
has to perform work on the particle, which, in turn, is 
partly dissipated in the fluid as heat, which has to be re- 
moved, and partly stored as potential energy in the har- 
monic potential. Macroscopically, both the work done 
on the system as well as the heat removed from it have 
to be positive in the NESS, while deviations from this 
macroscopic behavior arise due to fluctuations. 

The properties of the work fluctuations in this dragged 
particle model have been investigated experimentally and 
theoretically in the literature (see e.g. [J, 0, [!, H| ) , m par- 
ticular with respect to the validity of the so-called steady 
state fluctuation theorem (SSFT) (cf. @, E 0, S H ) ■ This 
theorem states that the probability distribution Tl T (p) 
of observing a particular value of a scaled dimensionless 
work value p over time r satisfies a certain symmetry 



relation, which can be formulated as [T(| 



n T (p) 

Hr(-P) 



(1) 



where = indicates the asymptotic behavior for large r 
and c is a constant. Eq. |T]) represents a refinement of 
the second law in that it quantifies the probability of 
observing temporary second law violations (negative p) 
in the NESS. The SSFT Eq. {TJ for the work fluctua- 
tions has indeed been validated in the dragged particle 
model if the noise is modeled as thermal Gaussian noise 
[HQ- Yet, if the noise is modeled as Levy noise the work 
fluctuations violate the SSFT and one can show that an 
anomalous fluctuation theorem [l3l[l4| holds. This result 
highlights the importance of explicitly verifying a general 
result, such as Eq. ([I]), in concrete models. 

In the present article we investigate the work fluctua- 
tion properties of the above mentioned NESS model when 
the noise from the environment has the characteristics of 
external asymmetric Poissonian shot noise (PSN), i.e., it 
is given as a sequence of one-sided Poisson distributed 
pulses with random amplitudes [lj|. PSN is a natural 
description of fluctuations in nature [TEI . [pol ] and is ubiq- 
uitous, e.g., in physics, electric engineering, and biology. 
Importantly, PSN generalizes the usual Gaussian noise, 
to which it converges in a certain limit, and is a paradigm 
of external noise: the dissipation and the noise originate 
from physically independent mechanisms 16]. This im- 
plies that the noise strength of PSN does not have to be 
related to the friction by a fluctuation-dissipation rela- 
tion, contrary to thermal noise in the case of Brownian 
motion. 

PSN allows us to investigate a number of important 
physical concepts in the context of nonequilibrium fluc- 
tuations: (i) Time scales: PSN introduces additional 
characteristic times in the system which crucially influ- 
ence its fluctuation behavior; (ii) Symmetry of the noise: 
PSN is usually asymmetric (one-sided); (iii) Singulari- 
ties: the asymmetry of the noise gives rise to an effective 
interaction between the noise and the potential, leading 
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to an effectively nonlinear system with singular features. 
Nevertheless, the work distribution can be calculated in 
analytical form in our model and is completely character- 
ized by its time scales. The interplay of these time scales 
leads to certain transitions in the behavior of the work 
fluctuations, as discussed in a recent Rapid Communica- 
tion [It]]. Here, we present the full theory of Ref. [l7| 
including inertial effects. Moreover, we consider the case 
where thermal Gaussian noise is superimposed on the 
PSN. This Gaussian noise models the effect of additional 
thermal noise on the particle, which might be relevant in 
an experimental realization of our model. 

The article is organized as follows. In the next sec- 
tion in we introduce the nonequilibrium particle model 
under the influence of one-sided PSN. Symmetries, time 
scales, and singularities of this model are discussed in 
Sec. IIIII The characteristic function of the work fluctu- 
ations including inertia is calculated in Sec. IIVI In the 
overdamped regime the characteristic function simplifies 
significantly and an analytical treatment is possible. This 
allows us to analytically determine the distribution of the 
work fluctuations using the method of steepest descent 
and investigate the fluctuation properties in Sec. |V] If 
the effect of inertia is non-negligible the work fluctuations 
exhibit time-oscillatory behavior, as studied in Sec. IVII 
Furthermore, in Sec. I VIII we discuss the combined ef- 
fect of PSN and thermal Gaussian noise on the work 
fluctuations. We close with some concluding remarks in 

Sec. Em] 



II. A NONEQUILIBRIUM STOCHASTIC 
PARTICLE MODEL WITH EXTERNAL PSN 

We consider a particle which moves in a time- 
dependent harmonic potential under the influence of fric- 
tion and external noise. The basic equation of motion for 
the position x(t) of the particle of mass m in the labora- 
tory frame reads (cf. 0]) 



mx(t) + ax(t) = -K(x(t) - vt) + 



(2) 



Here, the force —n(x(t) — vt) stems from a particle con- 
fining potential U(x,t) — k(x — vt) 2 /2 which is pulled 
with constant velocity v. For our purposes it is sufficient 
to restrict the discussion to one dimension. The parame- 
ter k denotes the strength of the potential, a the friction 
coefficient, and stochastic noise due to the environ- 
ment, to be defined below (cf. Eq. (|14[) ) . If the potential 
moves for a time period r, a certain amount of work is 
done on the particle, namely 



W T = -kv / (x(t) - vt)dt. 
Jo 



(3) 



In an experimental setup, where the harmonic potential 
can be induced e.g. via lasers 0, Eq. ^ is the mechan- 
ical work needed in order to move the potential in time 
r. If the potential is stationary (i.e. v = 0), no work is 



performed at all. This means that we ignore the work 
originating from the stochastic motion of the particle, 
which is attributed to the heat @, HH . 

The work is partly dissipated as heat due to particle 
friction, and partly stored as potential energy in the po- 
tential. In the steady state the mean value of the work 
(W T ) is positive, since we have to perform a macroscopic 
amount of work against the friction in order to maintain 
a NESS. This is basically a statement of the second law 
of thermodynamics (cf. jl9l I^V 

The mathematical treatment of our model Eq. @ is 
simplified if we transform to a coordinate system in a co- 
moving frame. Let us denote the position of the particle 
in the comoving frame by y(t) = x(t) — vt. The equation 
of motion for y(t) then reads 



my(t) + ay(t) = -ny{t) - av + £(t), 
and the work is given by 



W T = -kv I y{t)dt. 
lo 



(4) 



(5) 



Although the work is expressed in terms of the comoving 
coordinate y(t), Eq. © actually gives the work done in 
the laboratory frame [1, 0] , in whose fluctuations we are 
interested here. 

We consider the model Eq. under the influence of 
external PSN rather than the usual Gaussian noise, so 
that the motion of the particle differs in general from 
Brownian motion. PSN is specified by a sequence of delta 
shaped pulses with random amplitudes Tk and can be 
expressed in the form [l5| 



;(t) = J2r k 6(t-t k ), 



(6) 



fc=i 



where n t , the number of delta shaped pulses in time t, is 
determined by the Poisson counting process 



Pint) = y - J r e- xt . 



(7) 



The parameter A denotes the mean number of pulses per 
unit time (rate of pulses) so that there are Xt pulses oc- 
curring in the time interval [0, t] on average. In that in- 
terval, the time of the kth pulse is uniformly distributed. 
When a pulse occurs, its amplitude Tfc is sampled ran- 
domly from a distribution p(T). For p(T) we choose an 
exponential distribution 



p(r) 



-r/r 



To 



(8) 



where all amplitudes T are assumed here to be positive, 
i.e., the PSN that we consider is one-sided. 

The noise z(t) specified _according to Eq. © has the 
characteristic functional UM 



G z( t)[g{t)\ 



(9) 
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for a general test function g(t). If p(T) is given by the 
exponential distribution Eq. (jHJ) one obtains 



G z(t) [g{t)] = e X Io°°(T^m l ) 



di 



(10) 



This characteristic functional implies delta-correlated cu- 
mulants [2l[, which arise due to the delta shape of the 
stochastic pulses in Eq. ©: 



Cn (^1 j •••) tn) 



1 



n!Ar^(ti 



■ t 2 ) ■ ■ ■ S(t n -! 



\nG z{t) [g{t)] 
*»)■ 



g(t)=0 
(11) 



The first two cumulants, the mean and the variance of 
z(t), are therefore given by 



(z(t)) = xr , 



(12) 



(z(h)z(t 2 )) - (z(h)) (z(t 2 )) = 2Xr 2 6(h-t 2 ).{13) 

Here, the brackets (...) denote the usual ensemble av- 
erage, which, more precisely, represents a path-integral 
average with respect to the probability 7- > [z(i)] of a noise 
trajectory z(t) (cf. Q). Since it is convenient to have a 
noise with zero mean, we take in Eq. |(3J) for £(£) 



at) = z(t) - Ar . 



(14) 



This means that the noise in our model is considered to 
consist of a random shot noise part and a deterministic 
part that is equivalent to a constant negative drift force 
on the particle. Since the shot noise acts only one-sided 
by our convention, the noise £(i) is strongly asymmetric, 
even though its mean value is zero by construction. 

In contrast to the case of a Brownian particle, the 
noise is here external, which implies that the friction 
coefficient a is not related to the noise strength by a 
fluctuation-dissipation relation. In the absence of driv- 
ing (i.e. v = 0) the stationary distribution of the particle 
position is therefore not a thermal equilibrium distribu- 
tion in general. However, in a certain limit the £(i) of 
Eq. (TTJ} does become Gaussian noise, namely if we take 
the limits 



A — > oo, 



r ->o, 



while keeping 



Ar o = const. 



(15) 



(16) 



If this Gaussian noise is considered to be thermal, 
i.e., originating from an equilibrium heat bath, the 
fluctuation-dissipation theorem requires that this con- 
stant is given by 



(17) 



where (3 can be interpreted as the inverse temperature of 
the heat bath. This leads to 

(£(*)) = 0, (18) 
mat')) = 2a(3- 1 8{t-t'), (19) 



while all higher-order cumulants of are zero. In the 
limits of Eqs. (Tr5"|) — fTT)) . £(t) then represents the stan- 
dard Gaussian white noise of an equilibrium heat bath. 

The limit Eq. (fTS"]) means that both the waiting time 
between successive pulses (A -1 ) and their amplitudes 
(Tk) become very small, i.e., we have very frequent (at 
every 'time step'), independent, and very small pulses. 
This is the signature of Gaussian noise. The system can 
then indeed reach a state of thermal equilibrium if v — 0. 



III. SYMMETRIES, TIME SCALES, AND 
SINGULARITIES 



The system given by Eq. describes a damped har- 
monic oscillator with inertia, driven by a time-dependent 
force and by noise. For symmetric noise (e.g. Gaussian 
noise) the properties of the work fluctuations in this sys- 
tem are symmetric with respect to the direction of the 
pulling velocity v, i.e., the work distribution can only de- 
pend on the absolute value \v\. However, the presence of 
asymmetric shot noise, prescribed by Eqs. and |5]). 
breaks this symmetry. Since both the direction of v and 
the direction of the noise (given by the sign of r ) can in 
principle be either positive or negative, there are in total 
four different combinations of the two. It is easy to see 
that the case v > and To > is symmetric with the 
case v < and To < 0. Likewise, the case v > and 
To < is symmetric with the case v < and To > 0. 
Therefore it is sufficient to discuss only two of the four 
different cases. In the following we investigate the work 
fluctuations for v > and v < 0, while To > always. 

The oscillator itself is completely characterized by the 
two time scales 



r m = m/a, 



the inertial time, and 



a/ k, 



(20) 



(21) 



the relaxation time. Associated with the PSN are two 
additional characteristic time scales. Firstly, we have the 
time t\ defined as 



r\ = A 



(22) 



which is the mean waiting time between two successive 
pulses. Secondly, we can identify a time scale t p , defined 
as 



r„ 



(23) 



relating the mean amplitude of the pulses and the friction 
due to the driving, so that r p is the ratio of two indepen- 
dent external forces, due to the noise (To) and due to 
the driving (v), respectively. In total there are therefore 
four different time scales in our model: the inertial time 
r m , the relaxation time r r , the mean waiting time t\ and 
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t p . While the first two are intrinsic time scales of the 
oscillator, the latter two arise due to the particular type 
of external noise that we consider. We will see below 
that the four time scales r m , r r , t\, and t p fully specify 
the properties of the work fluctuations in our model, and 
that the interplay of these times induces transitions in 
the qualitative behavior of the work fluctuations. 

One familiar example of such a transitional behavior 
is due to the interplay of the inertial and the relaxation 
times, Eqs. (I2TJ1) and (j2"Tj) . respectively. Using r m and r r 
we can express the eigenvalues of the oscillator Eq. in 
the form 



^1.2 



4T m / T r 



(24) 



where the index 1 corresponds to the + sign and 2 to 
the — sign, respectively. We can thus formulate a critical 
condition (cf. Q) 



74 



(25) 



for the transition from real to complex eigenvalues. If 
Tm > T r /A the eigenfrequencies are complex and the in- 
fluence of the inertia on the dynamics manifests itself in 
time oscillatory behavior. Moreover, in a regime where 
T m *C T r inertia can be ignored and the system is ef- 
fectively overdamped. After neglecting the inertial term 
mx in Eq. the equation of motion in the overdamped 
regime can then be written as 



x(t) = - — (x(t)-vt) + -£(t), 
7> a 

or, respectively, in the comoving frame 



(26) 



y(t) = --y(t)-(v+^)+-z(t), (27) 
TV \ a J a 

where we have used Eq. (|14|) . We see that the subtracted 
mean value of Eq. (fT4")) has the effect of an additional 
drift force on the particle in Eq. (f27|) . Introducing an 
effective velocity v e defined as 



v e = v + Ar /a, 
allows us to write Eq. (|27|) as 



y(t) = —y(t)-v e + -z(t), 

T r a 



(28) 



(29) 



The one-sidedness of the shot noise z(t) leads to singu- 
lar features of the work fluctuations in the overdamped 
regime, as we discuss in more detail in the following. 



A. Singularities in the overdamped regime 



From the Langevin equation (|29p we can infer two 
important properties of the model in the overdamped 
regime. Firstly, upon averaging of Eq. (|29[) we obtain 



dt 



(!/(*)) 



1 



(y{t)) - v. 



(30) 



In the NESS the time derivative on the left hand side (lhs) 
is zero and the stationary mean position (y) is simply 
given by 



(y) = -vr r . 



(31) 



Then, using Eq. ([5]), we find the mean value of the work 
in the steady state 



(Wr) 



(y) dt 



av 2 r. 



(32) 



Clearly, (W T ) is always positive as required by the second 
law. These expressions for the mean position and mean 
work are the same as in the Gaussian case (see e.g. @) 
due to the zero mean of the noise in both cases. 

Secondly, we find that there exists a minimal value y* 
of the position coordinate. This can be understood if we 
consider Eq. (f2"9"|) without z(t), that is 



*(*) = — y(t) 



(33) 



In the NESS, where the lhs of Eq. (|33|) is zero, the particle 
will reach a position y* given by 



y 



(34) 



Since the shot noise part z(t) can only move the particle 
in the positive y-direction, this position y* is the min- 
imal position the particle can reach, i.e., y* represents 
a cut-off in position space. The origin of this cut-off is 
the asymmetric form of the noise Eq. (fl4")) . which is 
given as a superposition of a shot noise part z(t) acting 
in the positive direction only and of a drift part acting in 
the negative direction only. Without z(t) the particle will 
move until it reaches the position y* where the negative 
drift force is balanced by the positive spring force. Un- 
der the influence of z(t) the particle can then only reach 
positions to the right of y* (cf. Fig. |T|) . 

The presence of a position cut-off implies that the 
parabolic potential becomes effectively nonlinear: the ef- 
fect of y* is that of an infinite barrier in the potential. 
The effective potential for the particle is given by the har- 
monic potential U(y), truncated at y* (cf. Fig. []}. Since 
the effective nonlinearity arises due to the asymmetry 
of the PSN, our model exhibits an effective interaction 
between the noise and the potential. 

From Eq. (f5j) we find that the work rate w(t) is given 

by 



w(t) = -VKy(t), 



(35) 



i.e., proportional to the position. The position cut-off 
therefore implies a cut-off of the work rate w(t), so that 
we also obtain a cut-off value of the work in time r in 
the NESS, namely 



w: = 



y at = —vny t = 



(36) 
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U(y) 




y* <y> y 



FIG. 1: The truncated (effective) harmonic potential U(y) in 
the comoving frame for v > (regime (i)). The black bullet 
depicts the particle at its mean position (y) , while the effective 
infinite barrier for the particle position at y* is indicated with 
a dashed black line. For v > both the mean position (y) and 
the cut-off y* are always negative (cf. Eqs. (|31[) and (|34|l '). 



U(y) 


(a) U(y) 

v / 
< / 


k (b) 

| •<— / 
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y* o 


<y> y 


y* <y> y 



FIG. 2: The effective harmonic potential in the comoving 
frame for v < 0. Here, the mean position (y) is positive 
(Eq. (|31[) ) while the cut-off y* can be positive or negative 
(cf. Eq. pTlll. (a) Negative y* , if r p > r A (regime (ii)). (b) 
Positive y* , if t p < t\ (regime (hi)). 



We can also express the position and work cut-offs in 
terms of the time scales t\ and r p . Using Eqs (]2"2"]1 
and (|23p allows us to write y* of Eq. (|34[) in the form 




where a(v) denotes the sign function defined as 

( +1 , v > 
a(v) = I -1 , v < (38) 
[ , v = 0. 

The work cut-off of Eq. ([36]) then reads 

W? = <W r >(l + <7(«)Z£V (39) 

Due to T > the effective velocity, Eq. ([28]). obeys 
always v e > v so that y* < (y), i.e., the cut-off is always 
to the left of the mean position. Both work and position 
cut-offs have different characteristics depending on the 
sign of v and the ratio of r p and r r . We distinguish 
three different regimes which are important for the later 
discussion of the work fluctuations (Sec. fV|) . 

(i) v > (see Fig. Q}. In this case the mean position 
is (y) < and the cut-off y* of Eq. (|3"T|) is also always 
< 0. The work cut-off W* is then > (cf. Eq. J3S])) and 
denotes the maximal work done on the particle in time 
t. This is due to the fact that for v > the work rate 
w(t), Eq. (|33| . is larger for smaller positions y(t). 

(ii) v < and r v > t x (see Fig. HJi). Both y* and W* 
are then negative (cf. Eqs. (|3T|) and (|3^|) ). Moreover, 
the work cut-off is now the minimal work performed over 
time r, since for negative v the work rate w(t) of Eq. (|3"5|) 
is smaller for smaller positions. 

(hi) v < and r p < t\ (see Fig.[5jD). In this case both 
y* and W* are positive. As in (ii), W* is the minimal 
work performed over time r. However, a positive minimal 
work value W* implies that no negative work fluctuations 
can occur. 



Due to the existence of the cut-off values y* and W* 
the distributions of both position and work are non- 
Gaussian, unless one considers the Gaussian limit of the 
PSN Eqs. [fT5 ]) — ([XT ]) . In that limit y* -> -oo and 
W* — > ±oo. We remark that a similar work cut-off has 
been first observed in a Brownian particle model, where 
the moving potential is give n as a nonlinear potential of 
the Lennard- Jones-type [221 ] . 

The position and work cut-off are properties of the 
overdamped regime only. If the dynamics is influenced 
by inertia the cut-offs disappear, because the particle can 
actually 'overshoot' the minimal position y*, following a 
strong fluctuation of z(t) in the positive direction, due 
to its mass. However, this inertial effect is only relevant 
on short time scales and disappears if r ^> r m . In the 
asymptotic time regime therefore, in which we are mainly 
interested here, the fluctuation properties of the particle 
are described by the overdamped equation of motion and 
exhibit the singular features discussed above. 

IV. CALCULATION OF THE 
CHARACTERISTIC FUNCTION OF THE WORK 
FOR PSN 

An elegant method to calculate the probability distri- 
bution of W T has been presented in [lj, [H| • This calcu- 
lation is based on a theorem for a generalized Ornstein- 
Uhlcnbcck process 23], which is also applicable in the 
present case, if we write the equation of motion ^ in 
terms of a two-component system for the position y(t) 
and velocity u(t) of the particle in the comoving frame: 

m = «w, (40) 

«(*) = — -u(t) — y(t) - —v e + —z{t), (41) 

using Eqs. (|2"8"]) and (fTi]) . In vector-matrix notation this 
equation can be expressed as a linear first order differen- 
tial equation of the form 

y(t)=My(t)+ a + z(t), (42) 
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where we have defined the vectors 

y(t) - (v[t) 

z{t) = V frw j < 

while the matrix M. is defined as 



u(t) J ' 




M = 



1 

i i 



(43) 



(44) 



The derivation of the work distribution then proceeds 
along the following steps: 

1. The stochastic process y(t) of Eq. (|42p describes 
a generalized (two-component) Ornstein-Uhlenbeck pro- 
cess [l6j]. In order to determine the characteristic func- 
tional of y(t), defined by 



G y( t)[h(i)] = /ffltp j»y h(t)y 



(i)di 



(45) 



for a two-component test function h(t), we can then ap- 
ply the two-component version of the Caceres-Budini 
theorem [23| (see Appendix This theorem states 

that GyM [h(t)] follows from the characteristic noise func- 
tional G z ( t ) [k(i)] via: 



G y(t) [h(t)] = e iyoKo+Ja r k W d * Gz(t) [k(t)], (46) 



where k(i) is given by 



k(t) 



*a(*) 



3 AlT(s -* ) h(s)ds (47) 



and 



k = k(t = 0) = 



fc 2 (o) ; ' 



(48) 



The initial position yo and initial velocity uq of the parti- 
cle are contained in yo = (yo, uo) T - Due to the zero first 
component of z(t) (cf. Eq. ((43])) , the noise functional in 
Eq. gHD is given by G z(t) [k(i)] = G z{t) [k 2 (t)/m], where 



G z (t)[g(t)], the characteristic functional of the PSN z(t), 
has the exact form Eq. (|10[) . 

2. We obtain the characteristic function of the work 

G W M = {^) (49) 
by considering the particular test function (cf. [H, 03] ) 



( 1( /) = ( ^™e(r-t) 



(50) 



in the characteristic functional G y ( t )[h(t)], where 9(i) 
denotes the Heaviside step function. This can be seen 
immediately by substituting Eq. (|5"U|) for h(t) into the 
definition of G y ( t )[h(t)], Eq. (J45j), i.e. 

G y(i) [h(t)] =G Wr (q), (51) 
due to the expression for the work W T Eq. ([5]). 

3. Finally, the work distribution follows by performing 
the inverse Fourier transform of Gw T (q) ■ In the over- 
damped regime this can be done analytically using the 
method of steepest descent (see Sec. |V|) . 

We now evaluate the functional k(t) using Eq. (|47p . 
The matrix exponential exp{.M T } can be determined 
by diagonalization of the matrix M.. Using the partic- 
ular test function h(t) of Eq. (|5"0|) one obtains then the 
components of k(t) in a straightforward way 



kl(t) = i^9(r- 



v\ - v 2 

V\ 

ka(t) = — 9r- 

^1—2^2 ' V v l 

^2 



8(r-*)(^(l-e-(^)) 
(l-e^(--*))V (52) 



9(r-t) ^ (l-e^< T 



-*) 



(53) 



Substituting these expressions into Eq. (|4"6"| leads to an 
explicit expression for the characteristic function of the 
work: 



GwM 



exp \ iy k - 
1 



J VeK ; ( i (yiT + 1 - e^ T ) - -4 for + 1 - e"> T ) 

T m (z/1 - ^2) V^l ^2 



1 



1 At 



(54) 



where yo contains the initial conditions yo and uo and ()53[1 by setting t = 0. The work distribution now follows 
the components of ko are obtained from Eqs. (|52p and by performing the inverse Fourier-transform of Eq. (|54[) . 
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FIG. 3: The distribution of the particle position in the NESS, 
P(y) of Eq. (|57[) . One observes the divergence for y — » j/* 
if r r < ta. The cut-off values are given by yl = —1.25, 
2/J = —2.0, and w| — —2.5. Parameter values: r r = 1.0, 
v = 1.0, To = 0.5 [2|. 



which we have not been able to perform exactly. Nev- 



ertheless, Eq. (|54| can be readily used in order to de- 
termine the hierarchy of cumulants for the work, which 
reveal an oscillatory behavior if the influence of inertia is 
non-negligible. This is further studied in Sec. I VII In the 
next section we focus on the overdamped regime, where 
an analytical form of the work distribution can be ob- 
tained. 



V. ASYMPTOTIC FLUCTUATION 
PROPERTIES 

In the asymptotic regime r — > 00 inertial effects can 
be neglected and the work fluctuations behave as in the 
overdamped system. The characteristic function of the 
work in the overdamped regime is obtained by taking the 
limit m — > in Eq. (|5l)l . when the eigenvalues Eq. 
are given by v\ — r^ 1 and V2 ~ ~ T mi respectively. 
The characteristic function of the work Eq. (|5"1)) then 
becomes 



G wM ) = ex P \iy k + W (l - (l - e^) r r /r) + 1 ^ [ j^^^ - l) ^] , (55) 

where yo is the initial position and fco denotes the upper component of ko in the overdamped regime, which reads 
fco = —qav (l — e~ r / Tr ). The integral in Eq. ([55} can be expressed in closed form leading to an exact expression for 
the characteristic function of the work: 

Gw T (?) = (l + iqT v(l - e-^)) TTT ^ exp { lyo fc + iqW* T (l - (l - e^) 7) + ^ ( 1 + ^ " J ) } (56) 

I 



This expression simplifies if we consider the NESS of 
the system, where the initial position yo is drawn from 
the stationary nonequilibrium distribution. This distri- 
bution can be found by solving the Fokker-Planck equa- 
tion associated with the Langevin equation fZ§\i . We 
obtain then for the distribution of the particle position 
in the NESS (see Appendix |BT|) 

P(V) = Y^-^ [^(y ~ 2/*)) ^ 1 e-<*-* m W°, (57) 

where r(r r /r^) denotes the Gamma function [24| with 
argument T r /r\. One notices that the exponent in the 
prefactor of the exponential in Eq. (|57[) becomes negative 
if r r < t\, so that then P(y) diverges for y — > y*. This 
singular behavior is related physically to insufficient noise 
activation in the system when r r < t\. For, the time t\ 
is the average waiting time between two successive pulses 
of the shot noise. Therefore, if r r < t\ the system relaxes 
fast compared to the time scale of noise activation. In 
other words, the system relaxes 'too quickly' in between 
stochastic pulses and thus spends most of its time at the 



position that it would assume deterministically without 
the shot noise, which is y* . Consequently P(y) diverges 
for y — > y* (cf. Fig. [3]). 

We can now average the initial position yo in the ex- 
pression for the characteristic function, Eq. (|56p , over the 
stationary distribution Eq. (|57|) . The result reads 

where W* is given by Eq. (pjfj)) . The work distribution 
is obtained by carrying out an inverse Fourier transform 
of Eq. (|5"8|) and allows then an investigation of the work 
fluctuation properties in the asymptotic regime. 

A. Large deviation form of the work distribution 

The inverse Fourier transform of Eq. (|58[) can be calcu- 
lated analytically using the method of steepest descent. 
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It is convenient to introduce the scaled dimensionless 
value of the work p, defined by 



P ; 



(Wr) 



(59) 



The work cut-off W*, Eq. (p?9"|) . then gives rise to an 
extremal value p* of the scaled work p, defined by 



*- W * 1 , < s t p 



(60) 



This expression implies that the three regimes of the work 
fluctuations, discussed quantitatively in Sec. OH] (below 
Eq. l[3"9"])). correspond to three different regimes of p*, 
namely: (i) v > implies p* > 1. (ii) t> < and t p > t\ 
implies p* < 0. (iii) u < and t v < t\ implies < p* < 
1. 

The distribution function of p, denoted by U T (p), is 
obtained from the inverse Fourier-transform of Gw T , 
Eq. (58]) by 



Hr(p) 



G^e^^dc?. (61) 



2tt 



Using Eq. (|58|) on the rhs of Eq. (|6T|) as well as Eq. 
and Eq. (f36|) . we see that the distribution H T (p) can be 
written in the form 



2 /-oo 

Hr(p) = ^— / x(g)e Th <*>dg, 
where, the functions xil) an d h(q) are given by 
X (q) = (1 + iqr v)^( T +^~ 1 ) , 



(62) 



and 



h(q) = iqav 2 (p* — p) + 



1 



1 



- 1 



(64) 



respectively. 

For large r the integral in Eq. (|6"2"]l will be dominated 
by its saddle-point and can be approximately calculated 
using the method of steepest descent [26|]. Neglecting 
terms of order r -1 / 2 leads then to the saddle-point ap- 
proximation of IT T (p) in the form 



IMP) = 



27rV \h"(q)\ 



i8+Th{q) 



(65) 



where q denotes the appropriate saddle-point and 9 the 
angle between the deformed integration path and the real 
axis. Saddle-points are determined from the condition 
h'(q) = 0, which, with Eq. ijM]). reads 



(l + iqT v) 2 



0. 



avrx (p* - p) 

The solution of this quadratic equation is given by 



(66) 



Importantly, the square root y(p* 
ways real, since both p* 



p)/(p* — 1) is al- 
p and p* — 1 are either both 
positive (v > 0) or negative (v < 0), respectively. This 
means that the real part lZ(q±) = for all p and we can 
conclude that the appropriate saddle-point is 



>l- 



i 

T v 



1 - 



p* -p 



(68) 



since the original integration path (the real axis) can be 
deformed to go through q~- without crossing the pole at 
q = i/(T v). Furthermore 8 = in Eq. (f6"5)) . as is re- 
quired for a real probability distribution, since h(q) is 
real for purely imaginary q (cf. Eq. (|64|) ) and therefore 
the path of steepest descent through q- is parallel to the 
real axis (cf. the discussion in Sec. IV of |27|). 

Substituting the appropriate saddle-point Eq. into 
the expressions for x(q)j h(q), and h"(q) yields for the 
saddle-point approximation Eq. ([65]): 



n T (p) = 



i \A7n 



4vr \p* 



P -P 



T A P*-l 



exp 



r 



p -p 



p 



1 



(69) 



With Eq. (|60"f we see that the distribution II r (p) is com- 
pletely specified by the time scales r r , t\, t p , in addition 
to the measurement time t. 

One notices two different singularities appearing in 
Eq. (|6"9"]) . Firstly, the derivative of II T (p) diverges for 
(63) p ^ p* &s H'(p) oc \p* — p\~ x in leading order. This 
means that the approach of II T (p) to the cut-off p* has 
a vertical slope. Secondly, one notices that II T (p) itself 
diverges for p — > p* if x r /r A < 3/2. However, for large r 
this divergence occurs only in an isolated point and can 
be ignored. 

We find that Eq. (|6"9"|) yields an excellent approxima- 
tion of the distribution H T (p) at least for r > 10r r . This 
is shown in Fig. [?] where we compare Eq. (|69[) with a nu- 
merical inverse Fourier transform of Gw T and also with 
results from a direct simulation of the equation of mo- 
tion (|29|) using a Poissonian increment method [28| . For 
r = 5r r one notices a slight deviation in the negative 
tail for p < —1 between the saddle-point result and 
the numerical inverse Fourier transform. For r = 10r r 
the saddle-point approximation is in excellent agreement 
with the inverse Fourier transform over the whole range 
of p values. 

For very large r the distribution Tl T (p) exhibits the 
large deviation form [29( 



n T (p) 



,-rI(p) 



(70) 



with rate function 
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FIG. 4: Comparison of the analytic saddle-point approxima- 
tion for II T (p), Eq. (|69|l . with a numerical inverse Fourier 
transform of Gw T as weu as with results from a direct sim- 
ulation of the equation of motion ((29} . Parameter values: 
v = 1, Tr = 1, tx = 0.2, t p = 0.5. 



The rate function I(p) of Eq. (fTTj) is plotted separately 
for v > and v < in Figs. [5] and [6j respectively. In 
both cases one observes the strongly asymmetric shape 
of the rate function, which attains its minimum at the 
most likely work value, namely at p = 1, i.e., at the 
mean work value as expected. 

(i) v > (see Fig. [5} . The work cut-off is here p* > 1 
and denotes the maximal work done on the particle. As 
p — > p* the rate function I{p) approaches I(p*) with a 
vertical slope and ends at the finite value I(p*) = 1/t\ 
(cf. Eq. ([71])). The left side of I(p) is unbounded and 
becomes asymptotically linear for large negative p. From 
Eq. (|70)) it follows that the work distribution H T (p) de- 
cays exponentially for large negative p. 

(ii) v < and r p > t\ (see dashed and dotted curves in 
Fig. The work cut off is now p* < and denotes the 
minimal work done on the particle. I(p) approaches the 
end point I(j>*) = 1/t\ with a vertical slope as p — ► p*. 
Now, the right side of I(p) is unbounded and becomes 
asymptotically linear for large p, so that II T (p) decays 
exponentially for large positive p. 

(iii) v < and t p < r\. Here, the work cut off is 
< p* < 1 and no negative work fluctuations occur as 
becomes evident in Fig. [5] (solid curve). The behavior in 
the approach to the cut-off and for large p is as in (ii). 



B. Fluctuation theorem 

In the asymptotic regime r — > oo the work distribution 
Eq. (p9"} has the large deviation form Eq. ([70)) . In order 
to further discuss the fluctuation properties of work it 
is convenient to consider the dimensionless fluctuation 




FIG. 5: The rate function J(p), Eq. ([7I)l. for v > and r A = 1, 
plotted as a function of p for three different p* values. At 
every p* , I(p) assumes the finite value I(p*) = \/t\ = 1 here, 
which is approached with a vertical slope. For large negative 
p, I(p) becomes asymptotically linear. All curves correspond 
to regime (i). 



& 1.5 




FIG. 6: The rate function J(p) of Eq. J7TJ for v < and 
t\ — 1, plotted as a function of p for three different p* values. 
At every p* the rate function ends at the finite value I(p*) = 
1/t\ = 1 here, which is approached with a vertical slope. 
For large positive p, I(p) becomes asymptotically linear. The 
dashed and dotted curves correspond to regime (ii), where 
p* < 0. The solid curve corresponds to regime (iii), where 
< p* < 1 and no negative work fluctuations occur. 



functi 



fr( P ) 



n T ( P ) 
(w T ) ~n T (- P y 



■In 



where the constant a is defined by 

CtT\ 

a = ~FT> 
1 o 



(72) 



(73) 



and has the dimension of an inverse energy. In the limits 
Eqs. (fTS)) — (|T7|) . where the PSN becomes thermal Gaus- 
sian noise, the constant a is identical with the inverse 
temperature /3, so that f T (p) then agrees with the fluc- 
tuation function considered e.g. in [2|. 

The steady state fluctuation theorem (SSFT), Eq. (fTj, 
predicts that f(p) = lim r _ voo f T (p) = P- From Eq. f7D|) 
we obtain here instead 

f(p) = 2(p* - l)p 
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FIG. 7: The fluctuation function f(p), Eq. ([74]), plotted as 
a function of p for various p* values in regime (i). For p* > 
2, /(p) becomes negative if p £ [p + ,p*]. The approach to 
/(p*) is with a vertical slope, while /(p*) itself remains finite, 
Eq. (|75p (thin dashed curve). For large p*, /(p) approaches 
the straight line of the SSFT. 



FIG. 8: The fluctuation function /(p), Eq. ([74]), plotted as 
a function of p for various p* values in regime (ii). Using 
the antisymmetry of /(p) with respect to p we plot only p £ 
[0, — p*]. The approach to /(p*) is with a vertical slope, while 
/(p*) itself remains finite, Eq. (|75|) (thin dashed curve). For 
large p*, /(p) approaches the straight line of the SSFT. 



defined on the interval [— p*,p*] 30]. We see that the 
SSFT is violated in our model, even though we have 
identified a large deviation form of the distribution. As 
p — > p* the fluctuation function diverges like \p* — p]^ 1 / 2 , 
i.e., the cut-off is approached with a vertical slope (cf. 
Figs. [7] and [5]). However, f(p) itself remains finite at the 
cut-off and assumes the value 

f(p*) = 2{p* - l)p* - 2(p* - l^-^y- (75) 

We now characterize the behavior of the asymptotic 
fluctuation function f(p), Eq. (f?4"]) . for the three different 
regimes. 

(i) v > (see Fig. [7]). Here, f(p) of Eq. (J7J} has 
three zeros at p — and p± = ±2y/p* — 1 (due to the 
antisymmetry of f(p) we need not discuss the negative 
root p-). The zero p+ becomes significant when p* > 2, 
because then p* > p + and p can assume values in the 
interval [p + ,p*]. The crucial observation is then that f(p) 
from Eq. (|7i|) becomes negative for p € There 
exists therefore a parameter regime in which negative 
fluctuations of a certain magnitude are more likely to 
occur than corresponding positive ones (see Fig. [7J. In 
fact, since p* = 1 + t p /t\ we find that p* > 2 if t p > t\. 

This property is due to the strongly asymmetric tails 
of the work distribution H T (p) : the negative tail decays 
exponentially, while the positive tail decays more rapidly 
to the cut-off at p*, so that H T {p) < n r (— p) for p values 
in the vicinity of p* (cf . the discussion of the rate function 
in regime (i), below Eq. ([70]) ). It is important to note that 
despite the existence of a considerable negative regime of 
f(p), the second law is never violated: the mean value of 
the work is always positive, (W T ) = av 2 T. 



(ii) v < and t v > t\ (see Fig. [5]). In this case f(p) is 
zero only at p = and increases monotonically for p — > p* 
until the finite value f{p*) is reached. For large p* , f(p) 
approaches the straight line of the SSFT. 

(hi) v < and t p < T\. In this parameter regime 
no negative work fluctuations occur and the fluctuation 
function f T (p), Eq. (|T2|) can not be defined. 

For both v < and v > we observe a pronounced 
linear regime of f(p) for small p values. This is a gen- 
eral consequence of the large deviation form of Tl T (p): 
expanding the rate function I(p) around p = and tak- 
ing the ratio H T (p)/H T (—p) leads to a cancelation of the 
quadratic orders and therefore the linear term dominates 
up to order p 3 in the fluctuation function (cf. [29|). 



C. Gaussian limit and the SSFT 

In the limits Eqs. |15|)— (H7J PSN goes over into Gaus- 
sian noise. In this limit we should therefore reproduce 
previous results for a dragged Brownian particle in a 
parabolic potential (see e.g. 0, [![). In particular, we 
should obtain a Gaussian work distribution and the work 
fluctuations should satisfy the SSFT, Eq. (1]). 

In order to demonstrate this we expand the character- 
istic function of work, Eq. ([55]) in powers of To. Retaining 
terms up to second order in Tq yields 

lnGV x (<7) ~ iqctv 2 T 

- q 2 ry _«-*/*•)) -(76) 

This means that in the Gaussian limit of PSN the 
work distribution is Gaussian with mean (W T ) = uv 2 t 1 
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Eq. ([32]) . and variance 

{W 2 T ) - (W T ) 2 = 2Ty (l- - |(1 - e"^)) . (77) 

For thermal Gaussian noise the fluctuation-dissipation 
theorem requires that r 2 , = r^a/? -1 (cf. Eq. (fT7|) 1. In 
the asymptotic regime we thus obtain from Eq. ([77]) 

_ ( Wt) ^ s 2/3- 1 ( Wt ) , (78) 

i.e., the variance is proportional to the mean. The mean 
Eq. (|32|) and the variance Eq. (|78|) determine the distribu- 
tion of the work W T . For the distribution of the rescaled 
work p = W T / (W T ) one then obtains the Gaussian 

n T (p) - ,/ZS e -/5(^>( P -D 2 /4 ) (79) 

V 47T 

in the asymptotic regime. The saddle-point approxima- 
tion Eq. (|69p yields the same result, if one expands with 
Eq. ([60|) in powers of ta/t p 

ypHf^-^i^-l), (80) 

since t\/t p — > in the Gaussian limit. The fluctuation 
function f T (p), defined by Eqs. (j?2")) and (|75|) . then sat- 
isfies lim r _»oo f T (p) — P- We therefore confirm that, in 
the limit where the PSN goes over into thermal Gaussian 
noise, the SSFT holds. 



so that inertial effects become significant when m > m* . 
For Tji > m* the eigenvalues v\ and v 2 of Eq. (|24|) are 
complex conjugates and can be written as 

1/1,2 = n±iuj, (82) 

where we define 

7 ^ = V^m/Tr - I- (83) 

ZT m ^7"m 

Complex eigenvalues lead to a time oscillatory behav- 
ior of the position coordinate y{t), which oscillates with 
frequency lo. Time oscillations are also manifest in the 
work distribution, as observed in Q for Gaussian noise. 
In [J] it has been shown that for large r the oscillation 
frequency of the fluctuation function is the same as that 
for the position coordinate (i.e. lo). Since our model 
is essentially the same damped oscillator as was investi- 
gated in [U, only driven by a different noise, we expect a 
similar behavior of the work fluctuations for m > m* . 

Even though we do not have an analytical expression 
for the work distribution including inertia, we can in- 
vestigate inertial effects for the work fluctuations from 
the characteristic function of the work Eq. (|54[) via the 
cumulants of the work distribution c„ (t) , defined as [HI ] 

lnG WT (q) = f]^c n (T). (84) 

n—l 



VI. INERTIAL EFFECTS FOR FINITE TIMES 

For finite r inertial effects can be significant if r m is 
of the order of r r or larger. The eigenvalues v\ and 
determined by Eq. (fM)) . become complex if r m > r r /4. 
This condition can likewise be expressed in terms of a 
critical mass m* determined by Q 



rn 



or 
Ik' 



(81) 



We thus obtain the cumulants by calculating the deriva- 
tives of In Gw T (q) 



»(r) 



1 d n 



In G Wt (q) 



(85) 



q=0 



which yields from Eq. (|54|) expressions for Ci(r), the 
mean, and c 2 {t), the variance of the work distribution: 



(W T ) 



V K 



yo ( - (1 - e»n - Yl (i _ + Uq f J_ (1 _ ^rj _ 1 (1 _ 

^2 V\ I \V\ V 2 . 



1 ^2 
2 



771 \V 

r 2 



<1T 2 ) - (IT r ) 2 = 2^ f ™ ) f f 1 (l - e-(-*)) - 1 (l - e «*(-«>)) dt, 



t a Vtti(^i - i/ 2 ) 
and the nth-order cumulant reads 



(86) 
(87) 
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respectively. We note that in Eq. (|86|) the initial position yo and initial velocity uo of the particle appear explicitly. 
We now focus on the behavior of the mean and the variance. Using Eq. (f8"2")l allows us to rewrite the mean work 
Eq. (|86p in terms of trigonometric functions 



fl 2 + UJ 2 



,,2 



.2 , , ,2 



r(/i 2 + uj z ) + 2/i + e^ T 



sin(wr) — 2/xcos(wr 

2 2 



T m T r {(l 2 + LO 2 ) 2 

Likewise, one can express the variance, Eq. (|87[) . in the form 



(w 2 ) - (w T y 



2T 2 



- u (l + e' tT ^— sin(wr) — cos(o;t) 
sin(wr) — 2/i cos(wr) 



1 



(89) 



A(lu 2 t + n)(p 2 + lu 2 ) - 3u( M 2 - 3cj 2 ) - — (p 2 + to 2 ) 2 {I - e 2 ^ T ) 



t a (u 2 +w 2 ) 3 \2r m r r w 
+8we MT ((u 2 - w 2 ) sin(wr) - 2fiuj cos(wt)) 
-e 2/JT (^i(^/ 2 - 3w 2 ) cos(2wr) + w(3u 2 - uj 2 ) sin(2wr))] . 



(90) 



Both mean and variance show an oscillatory decaying 
behavior. The mean oscillates with frequency to while 
the variance exhibits oscillations with a superposition of 
frequencies u> and 2oj. For very large r we obtain from 
Eqs. © and (50]): 



A,(T) 



(a) 



A,(T) 



(b) 



(W T ) 

(W 2 ) - (W T f = 2T^t/t x , 



av 2 r, 

^2„,2. 



(91) 
(92) 



i.e., we recover the mean and the variance of the over- 
damped work distribution (cf. Eq. (f3"2"| and Eq. ([77]) ') 
respectively, as expected. The oscillatory behavior of the 
mean and the variance are shown in Fig. [5J where we plot 
the rescaled quantities 



Mr) 
Mr) 



(Wt) 
av z T 

(W 2 ) - (Wr) 

2T 2 v 2 t/t x 



(93) 
(94) 



using Eqs. (|9T[) and (|92[) . respectively. Both A\ and A 2 
converge to 1 in the limit r — > oo. 

Considering the higher order cumulants, one can 
also rewrite the nth-order cumulant Eq. in terms 
of trigonometric functions. Without performing the 
calculation, one notices that this would lead to de- 
caying oscillations with a superposition of frequencies 
uj, 2oj, 3oj, nu>. 



VII. SUPERPOSITION OF PSN AND 
THERMAL GAUSSIAN NOISE 

So far we have investigated the properties of the work 
fluctuation of the dragged particle of Sec. HH when the 
noise from the environment is given purely by external 
PSN. In this section we investigate the effect of addi- 
tional thermal Gaussian noise on the fluctuation proper- 
ties of the particle. The additional Gaussian noise takes 
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FIG. 9: Oscillatory behavior of (a) the rescaled mean Ai(t), 
Eq. (JS3J), and (b) the rescaled variance .4.2 (t), Eq. (IM)) . for 
the initial conditions yo = Uq = and various r m . In the 
limit t — > oo both Ai(r) and ^2(1") converge to 1. Parameter 
values: r r = 1, t\ = 0.2, To = 0.5. 



into account the effect of an additional equilibrium heat 
bath on the dynamics. This model can therefore rep- 
resent a dragged Brownian particle, which is subject to 
PSN. For a Brownian particle the thermal fluctuations 
and the friction have the same physical origin, namely 
the surrounding heat bath ('water'), so that the friction 
and the noise strength of the thermal noise are related 
via a fluctuation-dissipation relation. On the other hand, 
we assume that the PSN arises due to an external physi- 
cal mechanism that is independent of the thermal noise, 
which implies that the two types of noises are statistically 
independent. 

Gaussian noise is symmetric so that the particle will 
now be able to access all positions in the harmonic poten- 
tial. Both the position and the work cut-off of the purely 
PSN case, discussed in Sec. IIII1 are therefore expected 
to disappear. Our quantitative investigations start from 
the equation of motion for the comoving coordinate y in 
the overdamped regime 

y(t) = --y{t)-v e + -z{t) + -r){t), (95) 
r r a a 



which is Eq. 



with additional thermal Gaussian noise 
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xG v ( t ) [kit) /a] 



(99) 



(*?(*)) = o, 
(r,(t)r)(t')) = 2aP~ 1 5{t-t') 1 



(96) 
(97) 



where (3 is interpreted as inverse temperature of the equi- 
librium heat bath. Eq. (|97|) expresses the fluctuation- 
dissipation relation between friction and noise strength 
[32|. The characteristic noise functional of rjit) is given 
by [16] 



G v (t) [g(t)} = exp 



a 



g{tfdt 



for a test function git). The characteristic functional 
of y(t) can then be calculated analogous to the two- 
component case treated in Sec. IIVI using the theorem 
of Caceres-Budini [HI]. One obtains (cf. Eq. (j46]) ) [oH 



G yi t)[h(t)] 



^Woko — ive J k(i)dt 



G z{t) [k(t)/a] 



where k(t) is given by 



kit) 



,(ts)/T, 



/i(s)ds, 



(100) 



and 2/0 = vit = 0) as weu as = k(t = 0). The super- 
position of the two statistically independent noises zit) 
and rjit) in Eq. ([93]) leads therefore to the product of the 
corresponding noise functional in Eq. ([99]) . Substituting 



the noise functionals G z ( t ), Eq. (|10|). and G v m, Eq. ([98 
with the test function kit) /a as argument into Eq. 
yields 



{f°° 1 f°° 1 f°° ( 1 
iy k Q - w e / k(t)dt / fc(i) 2 di + — / . Fn , . . 

Jo a/3 J rxJo\l- i^k(t) 



1 dt 



(101) 



From this characteristic functional with kit) given by 
Eq. poop we can determine both the characteristic func- 
tion of the particle position and that of the work by 
choosing appropriate test functions hit). 



Distribution of the particle position in the 
NESS 



The characteristic function of the particle position is 
obtained from G y ( t ) [hit)] if we choose the test function 
hit) = hx5{t-tx) (cf. Appendix EH). The distribution 
of the particle position in the NESS is then obtained by 
carrying out an inverse Fourier-transform of the charac- 
teristic function. Details of this calculation are presented 
in Appendix IB 21 The result for the distribution can be 
written in the form of the convolution integral 



Pcpiv) 



P(y')P G (y-y'W, (102) 



where Piy) and Pdv) are the NESS distributions of the 
particle position in the purely PSN case and purely Gaus- 
sian case, respectively. Piy) is given by Eq. ([57]) and 
Paiy) is given by the Gaussian 



Pa(v) 




(103) 



The superposition of the two independent noises in the 
Langevin equation (|95p thus gives rise to a convolution 
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FIG. 10: The distributions of the particle position Pcpiy), 
Eq. (|102[) . and Piy), Eq. (|57|l for r r < t\ and various values 
of the noise ratio B. Piy) exhibits a divergence at the cut-off 
y* = —1.4, which is reduced to a maximum in Pgp(?/) due to 
the additional Gaussian noise. For B > the left tails of the 
various curves are Gaussian and the right tails exponential. 
Parameter values: t t = 1, t\ = 1.25, To = 0.5, v = 1, t t = 1. 



of the corresponding distributions. A closed form expres- 
sion of the integral Eq. (|102p is given in Appendix IB 21 in 
Eq. (IB1811 . 

In order to characterize the relative strength of the 
PSN to the Gaussian noise we introduce the dimension- 
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FIG. 11: The distributions of the particle position Pcp(y), 
Eq. (|102|1 . and P(y), Eq. (I57|l . for r r > t\ and various values 
of the noise ratio B. For small B values Pgp(j/) approaches 
P(y) (solid red curve), which exhibits a cut-off at y* — —2.0. 
For B > the left tails of the various curves are Gaussian 
and the right tails exponential. Parameter values: r r = 1, 
r A = 0.5, To = 0.5, v = 1. 



less noise ratio B, defined as 



B = 



(104) 



which represents the ratio of the noise strength of the 
Gaussian noise, 2a/ (3 (Eq. ([97])), and the noise strength 
of the PSN z(t), 2TI/t x (cf. Eq. fgfy). This means that 
large B indicates a dominant influence of the Gaussian 
noise and small B that of PSN. The limit case B = 
corresponds to the purely PSN case (see below). 

We plot Pcp{y) separately for r r < t\ and r r > t\, in 
the Figs . ITD1 and [TT1 respectively, for different values of B. 

In the case r r < t\ (see Fig. [TU)) the distribution P(y) 
of Eq. (|57|) exhibits a divergence at the position cut-off y* 
(cf. solid red curve in Fig. [10]). Under the superimposed 
Gaussian noise this divergence reduces to a maximum 
which shifts more and more to y* the smaller the values 
of B. At the same time one notices that in the approach 
to the maximum the curves of Pgp(jj) for B = 0.1 and 
B = 0.01 are basically on top of P(y) (cf. dashed black 
and dotted blue curves in Fig. [TO]) , i.e., the right tail of 
P GP (y) approaches P(y) for small B. The left tail on 
the other hand decays like a Gaussian 

In the case r r > t\ (see Fig. [TT]) there is no divergence 
in the distribution P(y). For small B, PGp(y) approaches 
the shape of P(y), yet without exhibiting a cut-off. This 
can be seen in the curve of Pcp(y) for B = 0.01 (cf. 
dashed black curve in Fig. [TT]) . which lies on top of the 
curve of P(y) (solid red curve) apart from a region in 
the vicinity of y*. For any B > the left tail of PGp(y) 
always extends beyond y* and decays like a Gaussian, 
indicating that the position cut-off vanishes due to the 
additional Gaussian noise. 

In both cases Pcp(y) becomes broader and broader for 
increasing B, i.e., stronger Gaussian noise (higher tem- 
perature). 



Work fluctuations 



From the characteristic functional G y (t) [h{t)], 
Eq. (jlOip . we obtain the characteristic function of the 
work Gw T (l) in a similar way as in the two-component 
case treated in Sec. IIVI by considering the test function 
(cf. Appendix |C"T|) 



h(t) = -qvnQ{r - t). 



(105) 



Substituting h(t) for h(t) into Eq. (|101[) leads to the char- 
acteristic function of the work 

GwM) = (1 + iqTov)^& exp {iqW? (l - — ) 



(W) 



3 T r 

27 



1 



(106) 



upon neglecting exponential terms in r and choosing the 
particular initial condition j/o = 0. In contrast to the 
purely PSN case, where the initial position yo has been 
sampled from the NESS distribution (cf. Sec. IIV|) . we 
consider here a fixed initial position for simplicity. In the 
r — > oo limit, in which we are interested in here, the par- 
ticular initial condition is irrelevant for the properties of 
the work fluctuations, because the work W T is extensive 
in r (cf. Eq. |5])). Therefore, a fixed initial condition 
yields the same result for the rate function as an average 
over the initial states. 

The distribution of the rescaled work p is then given by 
the inverse Fourier transform of Gw T (?) • This Fourier in- 
version can be performed as in Sec. [V] using the method 
of steepest descent. Details of this calculation are pre- 
sented in Appendix IC II The saddle-point q is now given 
by 



T Q v 

where r p is determined by solving the cubic equation 



(107) 



2Br p (r p ~l) 2 



P -P 
p* - 1 



(r p -l) 2 + l = 0, (108) 



where, p* denotes the work cut-off of the purely PSN 
model, Eq. ([60]) . Since there is no restriction on the pos- 
sible p values in the presence of the additional Gaussian 
noise, p* does not have the meaning of a cut-off here, 
but nevertheless appears as a parameter characterizing 
the PSN. An analysis of Eq. (|108|) shows that this cu- 
bic equation has a unique real root < 1 for all possible 
values of B, p*, and p (see Appendix IC 2|) . This par- 



ticular root (denoted by r p in the following) thus yields 
the correct saddle-point q, since the integration path in 
the saddle-point approximation can be deformed to go 
through q without crossing the pole at q = i/(T v) in 
the characteristic function Eq. (|106|) . 
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Using q of Eq. (|10T[) with the particular root r p , 
the saddle-point approximation of the work distribution 



J 



H GP (p) can be obtained in a straightforward way and 
reads (cf. Appendix IC ip 



ii gp (p) = 



x exp 



t/t x 



(1 



+ - 

1 9 



4tt(1 + B(l - r p ) 3 ) \p*-l\ 



P -P 



1 - 



- 1 



(109) 



Noting that p* = 1 + o{v)t p /t x , Eq. ([60]), U^ p (p) is 
completely specified by the times r, t Ti t\, t p and the pa- 
rameter _B. Furthermore, since r p < 1 always, the work 
distribution is real for all p values, so that there is no 
work cut-off as expected. 

We now investigate the behavior of the work distribu- 
tion in more detail in the asymptotic regime r — > oo, 
where H GP (p) assumes the large deviation form 



h^(p) 



(p) 



with rate function 



I GP (P) = ~ 



P -P 



r 2 B 

' J) 



(110) 



111) 



In Figs. IT2T— THlwe plot I GP (p) for various B values to- 
gether with the rate function I{p) of the purely PSN case, 
Eq. (|7Tj) . As before we distinguish three different regimes 
of the work fluctuations, namely (i) v > (Fig. [T^]); (ii) 
v < and t p > t\ (Fig. Q2]); (iii) v < and t p < t\ 
fFig. [T4|) . In all three cases I (p) is asymmetric around 
its minimum at p = 1 and becomes broader for increas- 
ing B. One observes that I GP (p) always extends beyond 
the work cut-off of the purely PSN case. This becomes 
particularly evident when comparing the B = 0.1 curve 
of I GP {p) (dashed blue curve in Figs. IT2T — Fl4|) with I(p) 
(solid red curve). I (p) lies virtually on top of I(p) on 
the unbounded side of I(p), but clearly deviates on the 
bounded side of I(p) and eventually increases monotoni- 
cally beyond the cut-off. Moreover, in case (iii) no nega- 
tive work fluctuations can occur in the purely PSN case 
due to the positive minimum work cut-off (cf. Eq. (|60p). 
but I GP {p) does indeed assume values for negative p. 
This means that in case (iii) negative work fluctuations 
arise due to the additional Gaussian noise. 

The behavior in the tails of I (p) for p — ► ±oo can 
be determined from the properties of the saddle-point 
q, Eq. (|107p . The saddle-point approximation namely 
implies that the slope of the rate function I GP (p) is pro- 
portional to the saddle point itself (cf. 29]), or more 
precisely (cf. Eq. (|C12|1 in Appendix IC ip 



tGP 



dp 



(p) = iqav 2 



-<r(i>)?, 



(112) 
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FIG. 12: The rate function I GP (p), Eq. (fTTTj) . plotted as a 
function of p for various B values together with I(p), Eq. (I71|l . 
in regime (i). The right tail of I(p) ends at the cut-off value 
p* — 2.5, where I(p*) = 1.0, while that of I GP {p) increases 
monotonically for any B > 0. Asymptotically the left tail of 
I GP (p) becomes linear and the right tail quadratic. Parame- 
ter values: t\ = 1.0. 



using Eqs. (fTU7|) and d23J). With Eq. (fTT2"j) one can de- 
termine the properties of the rate function as p — > ±oo 
as follows. 

First, we consider the case v > 0. As discussed in 
Appendix IC 21 r p has the properties that r p — > 1 for 
p — > — oo and r p oc —p for p — > oo if v > 0. It then 
follows from Eq. (|112p that I (p) has an asymptotically 
linear left tail and an asymptotically quadratic right tail. 
Therefore, from this behavior of the rate function one 
can conclude that the left tail of the work distribution 
Ii GP (p) for large r is asymptotically exponential, while 
the right tail is asymptotically Gaussian, if v > 0. 

For v < the situation is reversed because here r p oc p 



for p — > — oo and r p — > 1 for p — > oo (see Appendix I C 2p . 
The left tail of I GP (p) is thus asymptotically quadratic 
and the right tail asymptotically linear. Consequently, 
the left tail of H GP (p) for large r is asymptotically Gaus- 
sian, while the right tail is asymptotically exponential, if 
v < 0. 

One also notices in Figs. [T2T- THl that the rate function 
I(j>) of the purely PSN case represents a discontinuous 
limit of the combined PSN and Gaussian case. For any 
arbitrarily small non-zero B, both left and right tails 
of I GP (p) increase monotonically, while for B = 0, i.e., 
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FIG. 13: The rate function I GP {p), Eq. (fTTTj) . plotted as a 
function of p for various B values together with I{p), Eq. (JTTJ , 
in regime (ii). The left tail of I(p) ends at the cut-off value 
p* = -1.0, where 7(p*) = 1.0, while the left tail of i GP (p) 
increases monotonically for any B > 0. Asymptotically the 
left tail of J GP (p) becomes quadratic and the right tail linear. 
Parameter values: t\ = 1.0. 
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FIG. 14: The rate function I GP (p), Eq. (fTTTj) . plotted as a 
function of p for various B values together with I(p), Eq. (ffljl . 
in regime (iii). The left tail of I(p) ends at the cut-off value 
p* — 0.5, where I(p*) = 1.0, so that no negative work fluc- 
tuations occur in the purely PSN case. On the other hand, 
the left tail of I GP (p) increases monotonically for any B > 
and becomes asymptotically quadratic, while the right tail 
becomes asymptotically linear. Parameter values: t\ — 1.0. 



when I GP (p) = I(p), the rate function ends at a finite 
point for p = p* . This singular limit is due to the cut-off 
singularity in the purely PSN case. 

In order to further characterize the fluctuation prop- 
erties we define the dimcnsionless fluctuation function 



1 



In 



n? p (p) 



b(w T ) '~il gp (- p )> 

where the constant b is defined as 



b = 



1 + A 



(113) 



(114) 



The reasoning behind this definition is that b — ► a 
(Eq. ([75])) for B -> and that b -> (3 for B -> oo, i.e., we 
recover the fluctuation functions of the purely PSN case 
and the Gaussian case respectively, in the corresponding 



limits of the noise ratio B. In the asymptotic time regime 
the fluctuation function f GP (p) = lim^oo f GP (p) is 
given as 



f^(p) 



(p*-l) 2 (l + B) 



P 



y(r_ p + r p ) 



1 



1 



1 - r p 1 - r_ p 



(115) 



which is obtained by substituting the large deviation 
form Eq. piOp with the rate function Eq. (jllip into 
Eq. (|113|) . The fluctuation function Eq. (|1 15[) is a func- 
tion of p, p* , and B only. We discuss the behavior of 
f (p) separately for the three different regimes men- 
tioned above. 

(i) v > (Fig. [TS"]) , For large B, i.e., strong thermal 
Gaussian noise relative to the PSN, f GP (p) becomes lin- 
ear with slope 1 in agreement with the SSFT. For small 
B values, i.e., when the Gaussian noise is weak, f (p) 
approaches the fluctuation function of the purely PSN 
case, f(p) of Eq. (fT4")) . The important observation is that 
f (p) is negative for p greater than a certain p , where 
Po denotes the value at which f GP (p) intersects the p- 
axis. However, in contrast to the purely PSN case where 
the fluctuation function f(p) is bounded by f(p*) (cf. 
Fig. EJ), here f (p) — > — oo as p — > oo. This behavior is 
due to the tails of the work distribution II^ p (p) in the 
asymptotic regime as discussed above: the right tail of 
H GP (p) is Gaussian and thus decays more rapidly than 
the exponential left tail. Consequently II^ p (p) becomes 
increasingly smaller than H GP (—p) for increasing p, so 
that j {p) is monotonically decreasing for large p. The 
negative regime of the fluctuation function is thus even 
more pronounced than in the purely PSN case (cf. solid 
red curve in Fig. [15]), where the possible work values are 
bounded by the cut-off p* . 

(ii) v < with t p > t x (Fig. [TJ). Here, f GP (p) is al- 
ways positive for p > 0. For large B values one recovers 
the SSFT as expected. For small B values f GP (p) ap- 
proaches the fluctuation function of the purely PSN case. 
For any B > 0, f GP (p) is monotonically increasing due 
to the different behavior in the tails of the work distri- 
bution H GP (p), which behave now oppositely to regime 
(i): the left tail of H GP (p) is now Gaussian and decays 
faster than the exponential right tail. 

(iii) v < with t p < t\ (Fig. [TT)) . In this parameter 
regime the fluctuation function f(p) is not defined be- 
cause there are no negative work fluctuations. However, 
under the influence of the Gaussian noise negative work 
fluctuations do arise, as discussed above (cf. Fig. and 
one can discuss the properties of f GP (p)- For large B, 
J (?) approaches the SSFT as in the other two cases. 
For small B, f GP {p) becomes steeper the smaller the B 
value and in fact diverges as B — > oo. This is consistent 
with a vanishing probability of observing negative work 
(cf. Eq. I113p . because the fluctuation function, defined 
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FIG. 15: The asymptotic fluctuation function f GP (p), 
Eq. (|115|l . plotted as a function of p for various B values 
together with f(p), Eq. (|74|l . in regime (i). The different be- 
havior in the tails of II^ P (p) leads to a pronounced negative 
regime of f GP (p), which decreases monotonically for large p, 
while f(p) (solid red curve) ends at a finite value when the 
cut-off p* — 2.5 is reached. For large B the SSFT is recovered 
(solid purple curve). 



FIG. 17: The asymptotic fluctuation function f (p), 
Eq. (|115p . plotted as a function of p for various B values in 
regime (iii). f GP (p) increases monotonically for large p. For 
small B values the slope of f GP (p) becomes singular, while 
for large B the SSFT is recovered (solid purple curve). Here, 
p* = 0.5. 




FIG. 16: The asymptotic fluctuation function f GP (p), 
Eq. (I115|l . plotted as a function of p for various B values 
together with f(p), Eq. ([74}, in regime (ii). f GP (p) increases 
monotonically for large p, while f(p) (solid red curve) ends at 
a finite value when the cut-off p* = 1.0 is reached. For large 
B the SSFT is recovered (solid purple curve). 



by Eq. (TT3|), diverges for n^ p (-p) -> 0. 



VIII. CONCLUDING REMARKS 

(1) As one of our main results we have shown that 
the work distribution in our two models — the dragged 
particle with purely PSN and with PSN plus additional 
(thermal) Gaussian noise — exhibits a large deviation 
form but that the SSFT does not hold for general val- 
ues of the parameters. This result differs from previous 
studies of the noncquilibrium particle model considered 
here, where the noise statistics was taken as white Gaus- 
sian noise Q or white Levy noise (HI, respectively. 
In the Gaussian case the work distribution has a large 
deviation form, which satisfies the SSFT, while in the 
Levy case there is no large deviation form and the SSFT 



does not hold. The PSN case is in this sense intermediate 
between the Gaussian and the Levy cases. 

(2) One of the striking signatures of the work fluctua- 
tions in our two models is a parameter regime where the 
fluctuation function is strongly negative. This feature is 
due to the asymmetric tails of the work distribution for 
v > 0: if the noise is purely PSN the left tail decays expo- 
nentially while the right tail decays faster than exponen- 
tial to the positive work cut-off, where H T (p*) = 0. If we 
superimpose Gaussian noise the cut-off vanishes and the 
positive tail of the work distribution decays like a Gaus- 
sian. In that case the negative regime of the fluctuation 
function is even more pronounced and f (p) decreases 
monotonically for large p (cf. Fig. I15p . The unusually 
large negative fluctuations in our model could be useful 
for applications (cf. |33|). 

(3) We have considered the work fluctuations in the 
asymptotic regime of the NESS. Other well-known results 
of nonequilibrium statistical mechanics are the transient 
fluctuation theorems (TFTs) [ll|, H3 and the nonequilib- 
rium work relation (Jarzysnki relation, JR) [35| . which 
could be checked in our model in the case of PSN plus 
thermal Gaussian noise. The TFTs and the JR only ap- 
ply if the initial condition of the work measurement is 
sampled from a thermal equilibrium state (the JR re- 
quires in addition that the system equilibrates at the end 
of the work measurement). 

It is then important to note that in our model the ther- 
mal equilibrium state is not achieved simply by setting 
v = 0. In fact, the distribution of particle positions at 
v = 0, Pgp(u) of Eq. (|102[) . is evidently different from the 
Boltzmann equilibrium distribution. In order to gener- 
ate an equilibrium initial state one would have to assume 
that initially the system is only coupled to a thermal heat 
bath and that the PSN is only 'turned on' after the po- 
tential begins to move. 

(4) Several variations of the PSN that we consider here 
are possible. For example one could assume a pulse shape 
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different from the delta-peaks of Eq. ([6]) or consider other 
forms of the distribution of amplitudes p(T). In particu- 
lar, one could consider symmetric PSN which acts double- 
sided, e.g., by choosing a Gaussian distribution for p(T) 
or by adding two PSN of the form of Eq. ([6]), one with 
strictly positive amplitudes T and one with strictly neg- 
ative amplitudes, respectively. One would expect that 
those features of our model that are based on the asym- 
metry of the noise would disappear under such symmetric 
double-sided PSN. This will be discussed in more detail 
in a forthcoming publication. 

(5) Our model with purely PSN represents an effec- 
tively nonlinear system due to the infinite barrier in the 
potential induced by the noise. A genuinely nonlinear 
system modulated by a periodic field has been inves- 
tigated in [39j and has also been shown to violate the 
SSFT. 

(6) Our results might be relevant in the context of an- 
alytic frameworks for studying gene expression. Similar 
models as ours have been considered in [36j, 137], l38j to 
model the stochastic dynamics of protein concentration 
in a cell. Here, the PSN represents random bursts in 
protein production, which occur with an exponentially 
distributed number of molecules. 

(7) PSN occurs quite naturally in electric circuits, 
where the discreteness of the electron charge causes time- 
dependent fluctuations of the electric current. Our the- 
ory could thus be realized in an experiment similar to 
the resistance-capacitor dipole of IH| . If the Brown- 
ian Johnson-Nyquist noise is sufficiently weak compared 
with the shot noise one might be able to observe the 
strongly negative regime of the fluctuation function (cf. 
Fig. 

(8) We note that our theory could also be adapted to an 
experiment similar to that of Mahadevan et al [42j , where 
a lubricated rod of a hydrogel sliding on a soft vibrating 
substrate is considered as a model for biomimetic ratch- 
eting motion. Instead of the purely oscillatory external 
vibrations of [42| one could induce external asymmetric 
PSN, as considered here, which could then lead to work 
fluctuations with similar features as those presented here. 

(9) Our dragged particle model under PSN and Gaus- 
sian noise could also be realized experimentally by a 
micron-sized colloidal particle in water confined by a laser 
trap (43|, similar to the setup in [l], 0, H^|. Here, the 
fluid environment is the origin of the thermal noise, as in 
the case of the standard Brownian motion. Additionally, 
stochastic pulses in the form of PSN can be imposed on 
the colloidal particle using modulations of the laser beam. 
In this setup the laser beam gives rise to both the confin- 
ing moving potential and the shot noise. It would then 
be interesting to see if the negative fluctuation function 
predicted by our theory could indeed be measured. 
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APPENDIX A: THE CACERES-BUDINI 
THEOREM 

Consider the generalized multi-component Ornstein- 
Uhlenbeck process 

x(t)=Bx(t) + C(*), (Al) 

where x(i) and are n-component column vectors and 
B is an n x n square matrix. Let us assume that the 
characteristic functional of the noise C,{t) is known, 

G at) mt)}^(e^ m<(mt y (A2) 

where • denotes a scalar product and k(i) is defined as 
an n-component column vector ('test function'). The 
Caceres-Budini theorem then states that the charac- 
teristic functional of the process x(i) is obtained from 
G^( t ) [k(i)] according to 

G x{t) [h(t)]=e^G at) Mt)}- (A3) 

where xo contains the initial conditions and the func- 
tional k(i) is related to h(t) via 

/oc 
e (s " t)£?r h(s)ds, (A4) 

and ko = k(£ = 0). The process x(i) is therefore com- 
pletely specified by G^[k(i)], in particular, all cumulants 
of x(t) are obtained by functional derivation of the rhs 
of Eq. (|A3|) with respect to the components of h(t). 

This theorem follows upon substitution of £(t) = x(i) — 
6x from Eq. (TAT]) into Eq. flM]): 

Gem] = ^r^)-(*w-M*))d^ (A5) 

Partial integration yields 

G c(t) [k(t)] = e -^o^fr(-Ht)-B T Ht))Mt)^ 

= g-ikoxo ^ e i/ TO M*)-x(t)dt\ 

= e-^G x{t) [h(t)} : (A6) 
where in the second line we have set 

h(t) = -k(i) - B T k{t). (A7) 
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Consequently, k(t) is given as solution of Eq. (|A7|) . which 
reads 



k(t) = e~ w k - e 



fZ3 J / „sB 



e se h(a)ds. (A8) 



Here, the initial condition ko has to be chosen such that 
lim t ^oo k(t) = 0. 



APPENDIX B: NESS DISTRIBUTIONS OF THE 
PARTICLE POSITION 

In this Appendix we determine the distribution of the 
particle position in the NESS for two cases: (i) The noise 
is given by PSN. (ii) The noise is given by a superpo- 
sition of PSN and thermal Gaussian noise. We apply 
two different calculation methods. In (i) we solve the 
Fokker-Planck equation corresponding to the Langevin 
equation for y(t). In (ii) we use the characteristic func- 
tional G y m[h(t)]. 



1. NESS distribution for PSN 

The Fokker-Planck equation for the distribution p(y,t) 
corresponding to the Langevin equation (|29|) reads [la ]: 



Ft P(y,t) 



+X p(T)\p(y-T/a,t)-p(y,t)]dr, 

J —OQ 

(Bl) 

where p(T) is the distribution of the pulse amplitudes T. 
If we write the shift — T/a in the second term with the 
help of the shift operator exp{— (T/a)d/dy} and perform 
the integral using the exponential distribution of ampli- 
tudes Eq. (jSJ), we obtain 



Wt P(y,t) 



-X 



r o 



i 



a dy \1 + ^-d/dy 



p(y,t). (B2) 



A solution of this differential equation can be found as 
follows. In the stationary state we have dp(y,t)/dt = 0. 
Furthermore, since y(t) is a stationary process in a con- 
fining potential, the probability current vanishes as well. 
The Fokker-Planck equation (|B2[) therefore simplifies to 



P(») = 0.(B3) 



l 9 + „.)p (s) -Ar 



We now define 

M 



1 + l -id/dy 



P(y), 



(B4) 



i.e., P(y) is obtained from f(y) via 



P(y) = i 



a dy 



(B5) 



In turn, the equation for f(y) is only of first order: 



-y + v e ) (l + — !■) f(y) - X—fiv) = 0. (B6) 
T r J \ a dy J a 



The function f(y) is thus given as 

_ fy y' + v-r r , / 

Jy (»'+»cT r )r /i. V 



f(y) tx e J y (y +"= T «-) 



(B7) 



where y denotes the lower integration limit. From 
Eq. ([B2I we then obtain a result for the distribution of 
y in the NESS: 



Ply) oc M0Tr e ~h C dV 

a{y-y*) 



, (B8) 



where y* is the minimal value of the position in the steady 
state: y* = —v e r r (seeEq. ([M]) ) and the lower integration 
limit y has to be chosen y > y* for the integral to be well 
defined. Setting y = —vr r then yields 



a(y - y*) 



At*-— 1 



-ya/To 



(B9) 



The normalization constant can be calculated in a 
straightforward way and leads to the final result for the 
stationary distribution 



P(y) 



1 a ( a 



r(Ar r ) T \ r o 



-(y-y*)) e 



Atv-1 



-(y-y*)a/To 



(BIO) 

where r(x) denotes the Gamma function [23]. This result 
is used in Sec.[V]in order to sample the initial condition in 
the steady state. Furthermore, we note that the Fourier 
transform of P(y) is given by 



T{P{y)}= 



iuy 



(Bll) 



where u> is the Fourier variable conjugated to y. 



2. Steady state distribution for PSN with 
additional thermal Gaussian noise 



The characteristic functional of y(t) is defined as 

G vW [h(t)] = (expUj y(t)h(t)dt\y (B12) 

From G y u) [h(t)] the characteristic function of the particle 
position is obtained if we choose the test function 



hi(t) = M(*-*i), 



(B13) 
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since substituting Eq. (|B13[) into Eq. (|B12|) yields 

G y( t)[hi{t)] = (e iftlH ^)) = Gy(hi,h), (B14) 

where the rhs is just the definition of the characteristic 
function of the position y. Substituting h\{t) in the ex- 
pression for G y (t)\h(t)\, Eq. (jlOip . and taking the t\ — ► oo 
limit yields 

G y (hx) = fl-t^&i^ * e^v'Sp*'. (B15) 

The NESS distribution Pap(y) is the inverse Fourier 
transform of G y {h\). Noting that G y {hi) can be written 
as a product of the Fourier transform of P(y), Eq. IjBlip . 



and a Gaussian, we can express the distribution Pcp{y) 
as the convolution integral 

/oo 
P(y')P G (y-y')dy', (B16) 
-oo 

where P{y) is given by Eq. (|B10j) . and Pg{v) is given by 
the Gaussian 

pg{v)= vW e "^ 2 - (Bi?) 

The convolution integral in Eq. (IB16|) can be evaluated 
in closed form [4(| and yields: 



P GP (y) 



1 a 



r(Ar r ) r V 2^2 V 2a 2 



2X -(Ar r + l)/2 



e 2 



/3«r§ / Ar r \ /Ai> 1 1 



^k(2/ - y* 



/^(v-lrt-lW^hF, 



FoV^k 
' Ar r + 1 3 1 



2' 2 



(3n{y-y*) 



, (B18) 



where i_Fi denotes the confluent hypergeometric function 
of the first kind [13] and r(x) the Gamma function. 



APPENDIX C: CALCULATION OF THE WORK 
DISTRIBUTION FOR PSN AND GAUSSIAN 
NOISE 

1. Saddle-point approximation and the rate 
function 

The characteristic function of the work Gw T (l) is de- 
fined as 



(CI) 



We can calculate Gw T (<?) by substituting the particu- 
lar test function hit) of Eq. (| 1 05|) into the characteristic 
functional of y(t), defined by Eq. (|B12p 

G y(t) [h{t)] = ^ e -W>C*><^ = G W M)- (C2) 

Here, the last step follows due to definition of the work 
W T: Eq. ©. 

The work distribution U^ p (p), is then obtained as the 
inverse Fourier-transform of Gw T , Eq. (|106[) . i.e. 

^ Gw T (q)e- lqp{W ^dq. (C3) 



n? p ( P ) 



2tt 



After substitution of Eq. (|106[) in Eq. (|C3[) we see that 
n r (p) of Eq. (|C3|) can be written in the form 



n^(p) 



2 /• oo 
(XV T ' 



2tt 



X (g) e r/l(9) dg, 



(C4) 



where the functions xW) an( i h(o) are given by 
X (q) = (l + iqT v)^^ 

x exp | — iqp*av 2 r r + — av 2 T r ^j , (C5) 

and 

%) = iqav 2 [p* - p) - I 2 — 



If 1 



-1 , 



(C6) 



respectively. For large r the integral in Eq. (1C4|) will be 
dominated by its saddle-point q, which is determined by 
the condition h'(q) — 0. Straightforward algebra then 
yields 



Tnv 



(C7) 



where r p is determined by solving the cubic equation 
2Br p (r p ~l) 2 -f^-tf\(r p -lf + l = 0. (C8) 
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This cubic equation is analyzed in more detail in Ap- 
pendix EH 

The saddle-point approximation of II^ P (p) is given by 



n T GP b) 



av 



2^v wm 



x(q)e 



i6+Th(q) 



(C9) 



where 9 denotes the angle between the deformed integra- 
tion path and the real axis. The result for H GP (p) after 
substitution of the appropriate saddle-point q, Eq. (|C7|) . 
is presented in Eq. (|109p . 

From Eq. (|C9[) one can directly derive an expression 
for the rate function I GP (p). A comparison of Eq. (|C9[) 
with the large deviation form Eq. (1110|) yields 



(CIO) 



where h(q) of Eq. (|C6[) is more precisely given as = 
h(q(p),p), i.e., /i(q) depends on the dimensionless work 
value p via the saddle-point q(p) and via p directly. It 
then follows from Eq. (jClOj) that 

±i° p <p) - -|m?(p),p) 

= -JLh(q(p),pW(p) - ^h(q(p), tfpll) 

The first term vanishes due to the property of the saddle- 
point h(q(p),p) = 0. The slope of the rate function is 
therefore given by 



= iqav 2 , 



(C12) 



by differentiation of Eq. (|C6|) . The fact that the slope of 
the rate function is given by the saddle-point is a general 
result of the theory of large deviations (cf. (29]]). 



2. Analysis of the cubic roots 

The cubic equation (|C8|) can be written in the normal 
form 

^2+Xy+fl + lV+Iz^O, (C13) 



2BJ 



B) 2B 



where 



V -p 

p*-l' 



(C14) 



Although we could use Cardano's formula to investigate 
the behavior of the roots as functions of B and 7, we 
obtain the same information using simple calculus. We 
define the function 



/M^-( 2 + ^ + ( 1 + 2) r + 



1 



7 



2B 



(C15) 

Setting the derivative of f(r) to zero leads to simple ex- 
pressions for the location of the extrema of f(r), given 
by 



n = i, 



r2= 3 



Jb ^ 



Since the coefficient in the cubic term of f(r) is > 0, 
f(r) is monotonically decreasing for r — > — 00 and mono- 
tonically increasing for r — ► 00. This implies that, if 
T2 > ri, f(r) has a maximum at r\ and a minimum at 
r 2 . On the other hand, if r 2 < r l5 /(r) has a maximum 
at r 2 and a minimum at T\ . By substitution we find that 
f(r±) = 1/(2B), i.e., /(ri) is positive for all values of £>. 
Consequently, also /(r 2 ) > 1/(2B) if r 2 < ri. From this 
information about the extrema of f(r) we can conclude 
that f(r) must intersect with the r-axis either when ap- 
proaching the maximum at r\ = 1, if r-i > r\, or when 
approaching the maximum at r 2 , if r 2 < r\. Therefore, 
there exists always a real root f < 1 of the cubic equation 
(|C8|) . which is the appropriate root for the saddle-point 
q, Eq. (|107p . because the integration path in Eq. (|C4[) 
can be deformed to go through q without crossing the 
pole in the characteristic function Gw T , Eq. (|106p , at 
q = i/{T v). 

Moreover, the behavior of this root under a change of 7 
can be assessed qualitatively from this analysis. For large 
7, f remains in the vicinity of r\, so that f — ► 1 for 7 — > 
00. On the other hand, for negative 7, f remains in the 
vicinity of r 2 , which is proportional to 7 (cf. Eq. (|C16|) ). 
so that qualitatively f oc 7 for 7 — > — 00. These results 
can also be obtained from Cardano's formula for the three 
roots of Eq. ([C8]). 



Since7is given by Eq. (|C14|) , where v* = 1+<j(v)t p /t\, 
Eq. ([SO]) , we thus find that 

(i) v > 0: f — > 1 for p — > —00 and f oc — p for p — > 00. 

(ii) v < 0: f — ► 1 for p — > 00 and f oc p for p — » —00. 

This behavior of the root f determines the properties 
of the tails of the rate function I GP {p) due to Eq. (|C12j) . 
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